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Abstract. In |8], G. David suggested a new type of global minimizer for the Mumford-Shah functional 
in R 3 , for which the singular sets belong to a three parameters family of sets (0 < Si, S2, S3 < n). We first 
derive necessary conditions satisfied by global minimizers of this family. Then we are led to study the 
first eigenvectors of the Laplace-Beltrami operator with Neumann boundary conditions on subdomains 
of S 2 with three reentrant corners. The necessary conditions are constraints on the eigenvalue and on 
the ratios between the singular coefficients of the associated eigenvector. We use numerical methods 
(Singular Functions Method and Moussaoui's extraction formula) to compute the eigenvalues and the 
singular coefficients. We conclude that there is no (51,52,53) for which the necessary conditions are 
satisfied and this shows that the hypothesis was wrong. 
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The Mumford-Shah functional was introduced in ^3] as a tool for image segmentation. Let £2 be a 
bounded open subset of R™ (the screen) and g be a bounded measurable function defined on Q, (rep- 
resenting the image). The functional concerns pairs (u,K) where K is a closed subset of f2 and u is a 
function belonging to the Sobolev space H 1 ^ \ K). It is defined by 



where TL n ~ 1 {K) is the Hausdorff measure of co-dimension 1 of K. Let (u, K) be a minimizing pair of J, 
(wich always exists P3E|). The third term of the functional forces u to be close to g while, due to the 
second term, u has slow variation on fl \ K. Since no regularity is assumed for u across the singular set 
K, we may hope that for such a minimizer K is the hyper-surface across which g has great variations, i.e 
: the hyper-surfaces delimiting the contours of the image. 

The main difficulty arising in the theoretical study of the minimizers is the regularity of the singular 
set. First, let us notice that we may remove from K a set of H"" 1 measure which is not useful. Indeed, 
if (u,K) is a minimizer, there exists a smallest closed set K\ C K such that u € H l (Q, \ K\). The 
pair (u,Ki) is called a reduced minimizer of the functional. In dimension n = 2, Mumford and Shah 
conjectured that if (u, K) is a reduced minimizer for J and £l is bounded and smooth, if is a finite union 
of C 1 arcs of curves, that may only meet by sets of three, at their ends, and with angles of 2ir/3. 

This conjecture still resists but there exist partial results. In particular, A. Bonnet [2] showed that in 
the case n — 2, every isolated connected component of K is a finite union of C 1 curves. 

The crucial point introduced by A. Bonnet was a blow-up process, which leads to the notion of global 
minimizer of the Mumford-Shah functional. One way to prove the Mumford-Shah conjecture would be to 
get a complete description of all the global minimizers and then, if the global minimizers turned out to 
be simple, go back to the minimizers of the functional in a domain. The second step would be realized by 
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proving that if a minimizer is closed to a global minimizer (which is true via blow-up) then its singular 
set is smooth. 

Here we are concerned with the case n = 3. In this context, one may conjecture that the singular set 
of minimizer is a finite union of C 1 surfaces intersecting each other on a finite number of C 1 curves. 
Let us first describe the blow-up technique. Let (u, K) be a reduced minimizer of J, let x € and let 

t > 0, we set 

n ai t := r x (n-a!), 

9x,t(v) '■= t~ 1/2 g(x + ty), VyEttx,t, 
K Xtt := t-\K-x), 

u x ,t(y) '■= t~ 1/2 u{x + ty), Vyeo X)t . 

Then (u x ,t> Kx,t) is a minimizer of the modified functional J x ,t in £l x ,t where 

J x , t (v,G) := H n -\G)+ [ \Vv\ 2 + t 2 ( \v-g x>t \ 2 . 

Now, let us take a sequence (tfc)fc | and set (uk,Kk) '■= (u x ,t k , Kx,t k ) to simplify the notations. Such a 
sequence is called a blow-up sequence of (u, K) at x. It turns out that up to extraction, the sequence of 
sets Kk converges to a closed subset of R 3 . On the other hand, since the factor i -1 / 2 tends to infinity 
when t tends to 0, the sequence u& may not converge to a function having finite values. To overcome 
this difficulty, we have to subtract from Uk a function which is constant on every connected component 
of R 3 \ . More precisely, we have 

Theorem 1.1. There exists a closed subset Koo C R™, a function Uoo € Lj oc (R n ) and for each connected 
component V o/R" \ K^, constants (Pv>k)k such that up to a subsequence, 

Kk — ► Krya locally for the HausdorfJ distance, Uk — flv,k — ► u oo in L\ oc (V). 

Moreover, the limit pair (uoo,^oo) is a reduced global minimizer of the Mumford-Shah functional in R" 
(see Detinition M.l)) . 

Definition 1.1. Let K be a closed subset of R" and let u S i / 1 oc (R"). The pair (u,K) is a global 
minimizer of the Mumford-Shah functional in R" if the following properties hold. 

• For every open ball B in R", 7i n_1 (if n B) < oo and f B * K |Vw| 2 < oo. 

• For every open ball B in R™ ; for every pair (v,L) which satisfies the property above and such 
that 

a) L\B = K\B, b) v\ Rn \ B = U|r»\b, 

c) if x,y £ R™ \ (B U K) belong to a same connected component of R™ \ L, then they are also in 

a same connected component of R™ \ K , 

then 

(1) H n -\KC\B)+ f \Vu\ 2 < H n -\LnB)+ f \\7v\ 2 . 

Jb\k Jb\l 

From now on, we fix n — 3. Let us list types of reduced global minimizers {u,K) that are already 
known. For the first four types, the function u is constant on each connected component of R 3 \ K. 

(i) K = 0. 

(ii) K is & plane. 

(iii) K is the union of three half planes sharing the same edge and making angles 27r/3 with each 
other. 
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(iv) K is the half cone spanned by the edges of a regular tetrahedron from its center. In this case 
R 3 \ K has four connected components, each one being delimited by three infinite triangular 
faces. 

(v) (Cracktips) K is a half plane. Choosing coordinates such that K = {{x, 0, z), x > 0, z £ R}, the 
function u is defined by 

\lr 9 

u(rcos9,rsm9,z) = e\ — cos-+C, Vr > 0, < 9 < 2ir, 

y 7r 2 

where C is a constant and e = ±1. 
If this list was complete, from Theorem ll.il every blow-up limit of a reduced minimizer should be one of 
the listed global minimizers. Let us now describe an example of |8. for which this situation seems to be 
wrong. The whole argument is heuristic and is far from a proof. Let R > and C > 0, the domain is 
the cylinder fl = {(x, y, z), x 2 + y 2 < R, —R < z < R (see Figure 1) and g(x, y, z) := go{x, y)ip(z) where 

( C for < 9 < 2tt/3, 
go (r cos 9,r sin 9) := I for 2tt/3 < 9 < 4tt/3, 

[ -C for 47r/3 < 9 < 2tt, 

and ip is a smooth cut-off function satisfying < (p < 1 and 




for z > 1, 
for z < — 1. 



Let us now consider a minimizer (u,K) of the functional J associated to f2 and g. Since g = for 
small z, one may think that, for z close to —R, Kf\ {z = z } = 0. On the other hand, for z close to R, 
since g — go, we may suppose that K n {z = z } is close to the union of the three segments across which 
g jumps, i.e: {(r cos 9, r sin 9, zo) ■ < r < R, 9 = 2fc7r/3, k = 0,1,2}. We may then expect that K is 
the union of three regular surfaces meeting on a curve {^{t) : < t < 1} satisfying 7(0) = (0, 0, R) and 
7(1) = (xo,yo,zo) witlLZp > —R. See Figure 2 below. 



"x 




Figure 1. Cylindrical domain Vl Figure 2. Expected shape for the 

singular set K 



For < t < 1, a blow-up around the point 7(t) would lead to a global minimizer whose singular set is 
the union of three half planes sharing the same edge. These global minimizers should be of type (iii) and 
it is not necessary to introduce a new type of global minimizers if we suppose that the angles between 
the half plane are 27r/3. The situation is different when we consider a limit blow-up at 7(1). At this 
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point we expect a global minimizer (u*,K*) whose singular set is the union of three plane sectors with 
a common edge and that make angles 2tt/3. More precisely, intersecting K* with the unitary sphere S 2 , 
we obtain a set of three arcs of big circles M; cS 2 ,i = 1, 2, 3. These arcs are vertical, start at the north 
pole where they make three angles of 27r/3. Denoting their lengths Si, i = 1,2, 3, we obtain: 

3 

(2) K* = K 5lMM := R+ x (J M t . 

i=l 

No global minimizer with this kind of singular set is known. And we may thing that the previous list of 
global minimizer was noy complete. The hypothesis of G. David is the following ([8J, sections 76 and 80): 

Hypothesis 1.1. There exists a new type of reduced global minimizers (u*,K*) where by translation and 
rotation invariance there exists < 61,62,63 < tt such that K* = Ks lt s 2 ,s 3 and where the function u* is 
homogeneous of degree 1/2, i.e: 

u*(x) = \x\ 1/2 E*(x/\x\). 

Moreover, this new class of global minimizers may be generated by one of them, using translation, rotation, 
multiplication by —1 and addition of a constant. 

With this new type of global minimizers, the list of reduced global minimizers is closed. 

Remark 1.1. For n — 3, the homogeneity 1/2 is the natural homogenity for a global minimizer. In 
fact, if we suppose that u* is homogeneous of degree a then the equilibrium between the surface term and 
the Dirichlet energy term in leads to a = 1/2 or u locally constant. In our case, the homogeneity 
is impossible (we would remove any bounded piece of K and contradict 0},). ^° 9 e ^ a = 1/2; consider 
L = K \ B(Q,r) U dB(0, r) and v = in 5(0, r), then let R go to in (QJ to obtain a > 1/2 and let R 
go to +00 to obtain the second inequality. 

The paper is organized as follows. In section we set the notations. In section we find some 
necessary conditions satisfied by (Si, 62, 63) and £ if the hypothesis were true. SectionQJis devoted to the 
description of the numerical methods we have used to check these conditions. The numerical results are 
presented in section 

Acknowledgment. The author thanks Guy David for having proposed this work and for helpful support. 
For every fact concerning the Mumford-Shah functional, we refer to his book jH]. We are also indebted 
to Patrick Ciarlet and Monique Dauge for helpful informations on the Singular Functions Method and 
other related methods. 

2. Notations 

Let 6 = (61,62,63) be in (0,7r) , and let M±, M2, M$ C S 2 be three arcs of great circles starting from 
the north pole with relative angles 27r/3 and with respective lengths 5i, 62, 63. Without loss of generality, 
we will assume that 

Mi C Cx := S 2 n {(x, y,z) : y < 0, x = 0} , 
A/ 2 cC 2 :=S 2 n{(a;,y,z) : y > 0, x = v^y} , 

M 3 CC 2 := S 2 n : y > 0, x = - V3y} . 

Let P be the plane {(x, y, z) £ R 3 : x — 0}, in particular C\ C P. The open subset of S 2 : S 2 \ U 3 =1 Cj 
has three connected subdomains Oo, fix- f^o is symmetric with respect to P (i.e: (0, 1, 0) G fio) and 
Qi = R l (Slo), i = —1,1, where R denotes the rotation of angle 27r/3 around the z-axis (with the usual 
orientation) . 
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We will denote by Ss the domain 

3 

S S := S 2 \|jM 4 . 

i=i 

In the sequel, L 2 (S$), H 1 (S$) and H 2 (Ss) will denote the standard Sobolev spaces on Ss and A the 
Laplace-Beltrami operator on S 2 . We also define the closed subspaces : 

L 2 (Ss) := S^eL 2 (S s ) : J E = oj , 
We recall the following classical result: 

Theorem 2.1. Let S £ (0,7r) 3 . For any f in Lq(Ss), there exists a unique E in V 1 (Ss) such that 



(3) 



-AE = /, in Ss, 
<9„E = 0, on dS$. 



Equivalently, E is the unique solution in V (S$) of the variational problem: j VE • V5 = J /5 /or ewer?/ 
5 in ^(Sa). It is also the unique minimizer in V 1 (Ss) of the functional F(S) := 1/2 J IVS*! 2 — J fS. 
We will note E := — A^/ this solution. 

The operator A^ 1 is a compact symmetric operator on Lq(Ss). We will use spectral properties of such 
operators. In particular, Lq(S$) has an orthonormal basis of eigenvectors of A^. We will note fJ>i(5) > 
1^2(8) > • • • > the eigenvalues of — A^ 5 counting multiplicities and for k > 1, we set Xk(S) := l//j,k(S). 
Alternatively, these eigenvalues may be defined by: 



X k (5) := minmaxj/ |V5 2 | : S £ V k , [ S 2 = 1 



where the minimum is taken over all /c-dimensional subspace Vk of Lq(S$). 

When S — (0,0,0), it is well known that A/c(0,0,0) = 2 for k = 1,2,3 with associated eigenvectors 
(x, y, z) i— > x,y or z. In particular: 

(4) / |V5| 2 > 2 / S 2 VSe^(S 2 ). 



This property will be used at the end of section UH 

We will need some well known facts about the splitting in regular and singular parts of solutions to 
the Poisson problem with Neumann boundary conditions in a domain with corners. For this theory, we 
refer to 0, ^I] or ^21- Let us denote by & the end of M, for i = 1,2,3. The domain Ss possesses 3 
re-entrant corners of angles 2ir at £1, £2, £3 • 

Remark 2.1. In Pj], EU> only flat domains are considered. In order to prove Theorem \2.S\ we may use 
local smooth maps to transform the —A operator on Ss in an elliptic operator with smooth coefficients 
on a planar domain with a cut. In fact, it seems more natural to prove Theorem, \2.S\ directly. The main 
ingredients, Green formula, trace theorems, density results and use of polar coordinates do not change 
when one replace the planar domain with cuts by Ss. 

We begin in introducing a set of singular functions. 

Definition 2.1. For x € Ss, let ri(x) denote the geodesic distance on S 2 between the points £j and x. 
Using the usual orientation on S 2 , for x £ Ss in the neighborhood of £, 8i(x) denotes the angle at £j 
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between Mi and the smallest geodesic segments [&,x] (see Figure 3 below). We use (ri(x),9i(x)) as polar 
coordinates near to define 



( y/ri(xj\ ( 6i(x)\ 
Si(x) := 2 tan I I cos I — - — I ip(x/pi) 7 for x € Ss and i = 1,2, 3, 

where ip S C^°(R+,R) is a smooth cut-off function such that ip = 1 on [0, 1/2] and ip = on [1, +oo). 
The positive numbers pi,p2,p3 are chosen such that for {x € Mj : j ^ i} we have ri{x) > pi . In 
particular the functions Sj have disjoint supports. 




Figure 3. S A - 



For i — 1,2,3, the function defined above belongs to H 1 (S$). Moreover, this function satisfies 
homogeneous Neumann boundary conditions on dSs, we have J g = and Asi belongs to L 2 (Ss). In 
fact Asi = on {x : ri{x) < pi/2}. If S5 were a domain with a smooth boundary, then the quoted 
properties would imply: s, € H 2 (Ss). In fact, we have Si € H S (S$) if and only if s < 3/2. 

Theorem 2.2. Let f <G L 2 (Ss) such that j s f = 0, and /ei £ € i? 1 (S^) solves 

[{]A£ = /, in S s , 
d n T, = 0, on dSs- 

Then there exists £ 6 i/ 2 (Sa) and ai,a2,a!3 € R such that 

3 

S = £ + Qi,-Sj. 

i=i 

3. Necessary conditions 

Let (u*,^*), <5 and £* be as in Hvpothesis ll.il 
I s * condition 

Since (u*,K*) is a global minimizer, the function u* belongs to H 1 (B(0,r) \ K*) for every r > 0, thus 
£* G i? 1 (S5) and from Remark IT"T1 we have: 

(5) ^eff'fS^lO}. 
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2 nd condition 

Moreover, from Q with L = K*, we have for every r > and every v in H 1 (B(0,r) \ K*) such that 

v \dB(0.r) = u \dB(0.r)'- 

I |v«*| a < / |v<f. 

JB(0,r) JB(0,r) 

We deduce that u* is harmonic in R 3 \ K* and satisfies homogeneous Neumann boundary conditions on 
K. In term of £*, the last assertion reads 

/ — AE* = 3/4E* in S 5 , 

[ ' \ d„£* =0 on dS s . 

In particular 

(7) [ £* = 0. 
3 rd condition 

By the uniqueness assumption in Hvnothesis ll.il K* is unique up to rotation and translation. Thus, at 
least two of the lengths Si are equal. In the sequel, we will assume without loss of generality that S$ is 
symmetric with respect to P, i.e: 

(8) S 2 =S 3 . 
By uniqueness we also have 

(9) E*(x,y,z) = Y,*(-x,y,z), V(x,y,z)eS s , 

(10) or T,*(x,y,z) = -£*(-£, y, 2), V(x,y,z)eS s . 

4* n condition 

From ©,©,(0, we may apply Theorem O with / = -3/4E* and £ = E*. There exist E* <E H 2 (S S ) 
and a*, a£, a*. € R such that 

3 

(11) £* = E* + ]Ta* Sl . 

i=l 

Now let us return in R 3 , let 1 < i < 3 and let us define a blow-up sequence (uk 7 Kk)k of (u*, if*) at £j. 
Let be the line R£i and Pi be the half plane containing Mi and whose edge is Di. It is clear that the 
sequence (K^) converges to Pi locally for the Hausdorff distance. Let us study the Lj oc convergence of 
the sequence (iik)k- We denote by II the map R 3 \ {0} — > S 2 defined by n(ir) := sc/|a;|. Using the above 
decomposition, we have 

u k (y) = t l 7 1/ V(6 + i fe t/)-|e i A/c + ?yr /2 E*(n(6 + i fe y)). 

Let B an open ball ball of R 3 . Using the decomposition lUTf and the fact that for j : ^ i, Sj = in the 
neighborhood of £j, we have for t k small enough: 

Ufc(y) = \&/t k +y\ 1/2 (a* Sl (Tl(t; % +t k y)) + £*(Tl&+t k y))y 
= h(y)+ll k (y). 
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Now we introduce new polar coordinates in R 3 : (Ri(y), @i(y), Zi(y)) such that D< = {y : R L (y) = 0}, 
Pi = {y : @i(y) — 0} and the azimuth is uniquely defined by = 1 and Zi(0) — 0. We have 

6i(n(& + t k y)) = Qi(y) + 0(t k ), 
ri{ILfa+t k y)) = t k Ri(y) + 0(tl), 
\Zi/t k + y\ 1/2 = t- 1/2 + 0(l), 
uniformely in y € B. Thus, from the definition of Sj, we obtain 



I fe (y) = a*^R~{y) cos + G(Vh), 

For the second term, we use the fact that H 2 (Ss) is embedded in the Holder space C ' 7 for < 7 < 1, 
in particular choosing 7 > 1/2, we easily obtain 

n fc (y)-tfc 1/2 E*(e,) - o(l), 

— 1 /2 ~ 

uniformely in 5 and we deduce that Uk —t k converges to 

<(y) := a */^) C os^, inLL(R 3 )- 
We now use the following result [Hj 

Theorem 3.1. Every blow-up limit of a global minimizer is a global minimizer. 

The pair (u*, Pi) is thus a global minimizer and since we have supposed that the list of global minimizers 
was closed, the only possibility is that (u*, Pi) is a global minimizer of type (v). Consequently, we have 

\a*\ = VI < i < 3. 

The first consequence of this equality is to exclude the case symmetric (Eq. EJ. Indeed, in this case, 
we would have a* — 0. Thus S* is antisymmetric, this symmetry implies a* 2 — a\ and the additional 
information given by the last equality may be written: 

(12) K| = K|. 

5* n condition Let Lg A (Sa) and V^(Sa) be the subspaces of antisymmetric functions in Lq(S$) and 
V rl (S«5). If © and lTLfi|l are true, then, in particular 4/3 is an eigenvalue of the operator A^. Clearly, 
L 2 a (Ss) is stable by the operator — A^. Let us note — A^^ this restriction and n a .A{5) > ^2,a{5) > 
■ • ■ > its eigenvalues counting multiplicities. We set \ k ,A{5) ■= 1/ jJL k ,A{5)- For k > 1, we have 



(13) 



X k , A (S) := mmmaxj/ |VS 2 | : SeV k , f S 2 = ll , 



where the minimum is taken over all /c-dimensional subspaces V k of V^(S«5). 

The next results states that A 2lJ 4(<5) > 2. Consequently, ©, ©, and lHT)|l imply that 

(14) 3/4 = \ hA (S). 

Proposition 3.1. For 5 e (0,tt) 3 , A 2 ,a(<5) > 2 - 

Proof. Let 6 in (0,tt) 3 . We have V^(S (7ri7r)7r) ) D V 1 {S S ) (where S( OTi7ri7r) := S 2 \U 3 =1 C,). Thus, from fDty. 

A 2: a(<5) > A 2 ,A(7r,7r,7r), 

where for k > 1, 

A fe)A (7r,7r,7r) := min max/ / |V5 2 | : Sel4, / S 2 = 

v%cvi(s (ir)7riW) ), U S(5 Js 6 

dim Vfc — 
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We have \i,a{t, tt, t) = with associated eigenspace REi where Si = i on fij for — 1 < i < 1. It is not 
difficult to see that there is no other eigenvector in V^(S( 7ri7ri7r )) which is locally constant. 
Now, let E 2 7^ be an eigenvector associated to A2,a(tt, 7r, 7r). We split E 2 in E 2 = 5_i + So + Si, where 
suppiSj C fli. Let us fix i such that Si ^ 0. This function is a non constant eigenvector of — A restricted 
to fli satisfying Neumann boundary conditions, in particular J Q Si = 0. We set S :— Si off, so that 
supp S C fio- We also define S by 

S(x,y,z) = S(-x,y,z), V(i,!/,z) £ S (w) . 

We have to study two cases. 

case 1: S = S. Since S is symmetric, we can define E in F 1 (S 2 ) by 

)S(x,y,z) if (x, y, 2) G O , 

S(R(x,y,z)) if (x,y,z) eCL-i, 
S{BT l {x,y,z)) ii{x,y,z)€Oi. 

case 2: We set 5" := 5 — 5 ^ 0. This function is antisymmetric, in particular S' = on P. In this case, 
we set: 



Y,(x,y,z) 



S'(x,y,z) if (x,y,z) eQ , 

S'(-x',y',z') where (x',y',z') = R(x,y,z), if (x, y, z) e f2_i n {y < 0}, 

S'{-x l ,y , ,z') where = E-^y,*), if 2) G Qi n {y < 0}, 

otherwise. 



In both cases E G V rl (S 2 ) \ {0}, satisfies / g2 E = and 

/ |VE| 2 = A 2 ,A(7r,^,7r) / E 2 . 
Js 2 Js 2 

Thus, from we conclude that \2,a{k, tt, t) _ 2. □ 

Complete problem. We now collect the necessary conditions obtained in this section. If Hvnothesis ll.il 
is true, from 1)51617181101121 and ITijl ■ then there exists 5 = (Si, 8 2 , 82) G (0, ir) 3 such that: 

(15) 3/4 = Ai,a(<5) = min |y |VSf : 5 G F X (S 5 ), J S 2 = 1, 5 antisymmetric | . 
Moreover, letting 

E(<5) G argminjy \VS\ 2 : S e V 1 (S 5 ), J S 2 = 1, 5 antisymmetric j , 

then the singular coefficients ai(<5), a2(<5), as(<5) such that E(<5) — J2i<i<3 oa(8)si G i/ 2 (S,5) satisfy 

(16) M<5)| = |a 2 (5)|. 

In the sequel, we give numerical evidences showing that there is no pair (5, E(<$)) satisfying both ltl5|l 
and iflfijl . The conclusion is that Hvnothesis ll.il is false. 



4. Numerical methods 

The general method is the following. Let ft. > 0, we choose a subdivision = 8$ < 8^ < ■ ■ ■ < 8^ < 
= 7[ . j satisfying 8l +1 - 8% < h for < k < N. Then, for every 8 h := , 8^ ,8%J (1 < ki,k 2 < N), 
we compute numerical approximations of Xi ! a(8 h ) and of the coefficients cti(S h ). Finally, we use these 
values to test the validity of equalities ltl5|l and iflfijl . 
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We use a Galerkin method to approximate \i,a(5 )■ More precisely, we set 
(17) Af iA (5 h ) := minJ f \VS 2 \ : S G V 1 ' 





(18) Z h (5 h ) G argmin<j / |V5 2 | : S G V 

ls 5h 

where V h (5 h ) is a finite dimensional subspace of V 1 (S S h). This space is chosen great enough such that 
we may hope that A (S h ) and T, h (Sh) are close to Xi,A{S h ) and S(5h). Typically, F ,l (£' 1 ) is the space 
of P 1 finite elements constructed on a triangular mesh of S$h of size h. 

Remark 4.1. We use the same letter (h) to denote the step size of the subdivison Sq < ■ ■ ■ < 6%f +1 and 
the mesh size of the triangular mesh of Sgh . These sizes could be different but they are actually equal in 
the numerical computations below. 

Let (T h )h>o be a family of regular meshes of 5 2 composed of geodesic triangles and with mesh size 
h. We assume that the edges of T h do not cross the geodesic segments C\, C2, C3. We also assume that 
R(T h ) = T h and that T h is symmetric with respect to P. This last symmetry is imposed in order to 
work with antisymmetric functions. We choose the subdivision = Sq < 5% < ■ ■ ■ < 5% < 5% +1 = it such 
that (0,sini5f ,cos(5f)o<i<Ar+i are the coordinates of the nodes of belonging to C\. 
Let h > 0. From now on, S h — {8^,5% ,5% ) and to lighten notations, references to S h will be omitted. 
Let (pi)\<i<M h be the set of nodes of T h and let Vh be the polyhedral domain of vertices {pi)i (the 
boundary of the convex hull generated by (p»)i). Recall that II is the projection of R 3 \ {0} on S 2 . This 
map defines a bijection from Vh onto 52, let us note II -1 its inverse. 

Now let ((Pi)i<i<N h be the set of continuous functions defined on Vh \ u| =1 II _1 (Mfc) such that the 
restriction of (p\ on each face of Vh is linear and such that there exists 1 < < M h such that 
<fi(pj) = 1 for j = otherwise. 

Finally, for 1 < i < N h , we set $ :— (p 1 } o II" 1 and we define the space of P 1 finite elements on S$h to be 

W h :=span{^ : 1 < i < N h } . 



And then 

V h :=\<p h e W h 



: / ^ 
J s 5h 



Remark 4.2. In general, the elements ofW h are not continuous across the geodesic segments M\, M2, M3. 
Let us also stress that we have N h > M h . Indeed, if pj belongs to Mi and pj 7^ £1 then there exists 
ii 7^ 12, such that f^iPj) = <Pi 2 (Pj) = 1 and supp^ C supp</? l2 C Oi. 

Remark 4.3. The constant functions belong to W h (indeed, X^Li <Pi = ^) an & V h is the orthogonal of 
1 for the L 2 inner product. 

Let / G L 2 (S sh ) and S : — A^ 1 ^. Since S does not necessarily belong to H 2 (S(S l1 )), the classical 
convergence rates obtained for the approximation by P 1 finite elements for a similar problem on a smooth 
domain are not valid here. In fact, for quasi uniform meshes, there exists c > 0, such that 

s mi n ft IS"' 1 - S\ H i(s sh ) > cmax|ai|\/ft, 

where ai, 02,0:3 are the singular coefficients of S. This conclusion holds for S := £(#/j). To overcome 
this difficulty, we add the singular functions to the space of finite elements. Namely, we set: 

(19) V h := V h 8 span{s 4 : 1 < i < 3}. 
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This method is called Singular Functions Method (see 011312 f° r a review on such methods). The usual 
approximation rates (valid for smooth domains) are recovered. 

nun h \S h -E\ HHS5h) < Ch. 

A classical result (see ^H] , for example) concerning the approximation of the eigenvectors of an elliptic 
operator by Galerkin methods leads to 

\\U(6 h ) - \ h A(t h )\ < Ch 2 , 

\Z\S h )-Z(S h )\ HHSih) < Ch, 

\X h (6 h )-Z(S h )\ L * {Seh) < Ch 2 . 

For the approximation of the singular coefficients on (8 h ) , we use an extraction formula of Moussaoui 
(see also j4]). We first introduce the dual singular functions 

Definition 4.1. With the notations of Definition \2.1l we define 

1 Odx) 

SAx) := ; — cos — - — ib(xlpi), for x E S$ and i= 1,2,3. 

2tan v /^)/2 2 

Now, forl<i<3, let p % G V l {S s ) be the variational solution of 

-Ap l = ASi, in Ss, 



(20) 
Finally, we set 



dpi = 0, on dSg. 



Pi := Si + p t , for 1 < i < 3. 



Remark 4.4. We have ASi = on {x : ip(x/pi) = 0}, so AS* is smooth and pi is well defined. The 
function pi does not belong to H 1 (Ss) (we only have pi € H S (S$) for s < 1/2). 

Theorem 4.1 (Moussaoui, ^s3)- / ^ Lq(Ss), S := Ajf s f and let ol\, 012,0.3 be the singular coeffi- 
cients of S. Then 

cti = — ( p l (x)f(x)dx. 



In our case, the singular coefficients on{8 h ) are obtained by the formula above with / :— \\ A (S h )T,(S h ). 
In order to get numerical approximations of these coefficients, we first compute an approximation p\ € V h 
of the functions pi . We have 

\Pi-Pi\ L * < Ch 2 , l<i<3. 
Then, we set p\ := Si + p\ and finally: 

a>?(6 h ) := \ h (5 h )- f p^(S h ). 

The numerical convergence rate is given by 

\a';(S h )-a t {S h )\ < Ch 2 , l<i<3. 

Figure 4 represents the error e(h) on the computation of Ai j a(^' 1 ), cei(5 h ) and a.2(5 h ) for 1/193 < h < 
1/5 and S h = (tt/2, tt/2, tt/2). The "exact" solution is obtained with h = 1/320. 

For the choice of ip (Definitions 12.11 and 14.111 , it is sufficient to have a C 2 function (we have used 
a piecewise polynomial function). The main obstacle for the accuracy of the numerical computations 
turned out to be the restriction on pi ( Definition 12. 1|) . If 8 h = {81,82,82) is such that one of the 8% is 
close to or ir, then we have to choose a very small pi. Consequently, the function ASi has great values 
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and we need a fine mesh to get an accurate approximation of pi. For this reason, we have worked with 
this method for S h £ (0.1, 3. 04) 3 . For other values of S h , we use a method based only on finite elements 
(without singular functions) described below. 




Figure 4. Full line: |A^(5 ft ) - Xi, A {S h )\. Dashed line: \a^(8 h ) - ai(5' l )|. Dotted line: 
\c%(6 h ) - a 2 {S h )\. Xi,A(S h ) ~ 0.795, ai (5 h ) « 0.54 and a 2 (S h ) ra 0.27. 

Let us fix S € (0,7r) 3 . Let ^f(S) (resp. H~(S)) be the west (resp. east) trace function of S(<5' 1 ) on Mj. 
From the definition of aj(5), we have 

oi) . im Si™) _ ± „ i(j) . 

2 

(The sign ± depends on the orientation choice of Definition 12 .11 1 

Now let h > 0, such that it is possible to set 6 h — 8. Let 1 < i < 3, and £q, ■ ■ ■ ,S,^h be the nodes of 

the mesh Th belonging to Mj. We suppose that the third coordinate of the sequence (£ )& is decreasing 
(in particular £q = (0,0, 1) and = £j). We replace y' 1 by V h in 11811 to compute an approximation 
'S EF (8 h ) of £(<$) and we define a new approximation of the coefficient a^((5) inspired by 12111 . 

(22) <bf(^) := ; 



4 tan 



2 

where T,^ Fi (8) and T, FFi (S) are the traces of T, EF (S h ) on Mj. In fact, the coefficients ofy EF (5) do not 
converge to cti(8), when /i goes to 0. However, since we are concerned with the ratio |a§ ( ( 5' i )|/|a ; i'(<5 ft ')|, 
it turns out that the method makes sense. During numerical experiments, we have observed that, if we 
consider a family of quasi-uniform meshes 7^ such that, the family of rescaled meshes 1/h x (7/,. — £j) 
tends to a fixed mesh % of the plane {x £ R 3 : • x = 1}. Then 

(23) lima^W = 

where q is a constant depending on %. We did not prove this claim. 
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In our numerical study, the mesh has the same shape in the neighborhoods of the three points £1, £2, £3- 
Consequently, we have c\ = c 2 = C3. Thus we may consider that |a2._EF(^)l/l a i ef($)\ is an approximation 
of |«2 (<5) I /|cki (5) I . We have compared the numerical convergence of both methods for this ratio. Since 
the exact ratio is not known, we have used \a% (8)/afl (8)\ on a finer mesh (h « 10~ 3 ) to evaluate the 
error. The numerical convergence results for the two methods are given Figure 5. 




Figure 5. Upper curve: 



Wi,S,Ef\ 



|a2,5|/|Q!l,<J 



Lower curve: 



\ a l,8\ ~ \0t2,s\/\0tl,s\ 



. For 1/193 <h< 1/10, S 1 = tt/6 and S 2 = 5tt/6. 



The interesting fact is that they do converge to the same limit. We will see in the next section that 
the zone for which the equalities l(T5|l . iflfijl are the more close to be true is the neighbourhood of (Si = 0, 
8 2 = 7r). For this reason, the comparison have been realized for Si = n/6 and 82 — 5n/6. 



5. Numerical results 

We have computed T, h and a\ with the numerical method described in the previous section for h 
1/40. We obtain a curve of approximate solutions of l(T5ll : 8% — f h (Si) (see Figures 6,7). 

Remark 5.1. Let us note that for Si — ir andS2 = 0, the corresponding singular set is a half plane. In this 
case, a global minimizer do exist: the cracktip (type (v)) and it is natural that lim^ ) ^ia(^Ii^) = 3/4. 
Actually, for (81,82) close to (ir,0), we observe that S' l (<5) is close to the trace of the cracktip on S 2 . 

We also have exact values for (^1,^2) = (0,0) and (81,82) — (0,7r) for which Ai,a(0, 0) = 2 and 
Ai,a(0,7t) = 21/16 = 1.3125. In the first case, the corresponding eigenvectors are H{(x,y, z) 1— ► y}. In 
the second case, §5 has two connected component, the space of eigenvectors associated to 21/16 is RS 
where S = on the small connected component and S(x, y, z) :— (cos^) 3 / 4 sin(3/46*) on the big connected 
component. (The spherical coordinates (<f>,6) are defined by \0\ < 2ir/3, \<f>\ < it/2 and (x,y,z) = 
(— cos (j) cos 9, — cos (f> sin 0, sin <fr) .) 

In order to check the condition l|T6jl . we compute the approximate coefficients a'KS' 1 ) for h « 1/40 and 
0.14 < 81,82 < 3. To complete the study we have computed the alternative approximations of the ratio 
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Figure 6. Values of \ h s for h = Figure 7. Curve tff := f h (6' 2 l ), h = 1/160. 

1/40 and 0.14 <Si,S 2 <3 . 



KjI/KjI given by K )BF (<J h )|/|af )BF (<J' , )| for ft « 1/160 on the curve : < $ < tt} 

(Figure 8). In both cases, we have 

\a%(5)\/\at(6)\<Q.8, 

for any couple (61,62) of the discretizations. This inequality contradicts llllit . 

For both methods, the numerical error is less than 1/100 (see Figure 5). This numerical error is small 
compared to the distance between 0.8 and 1. We conclude that there is no value < 61,62 < tt for 
which an d QHt are both satisfied. Consequently, we are convinced that Hypothesis 11.11 was wrong. 

Remark 5.2. We observe that a>i(5) converges to when 61 tends to so that solutions of l|16H do exist, 
but in this case 61 < /(J2) (ind \i,a(6) > 3/4. 

Remark 5.3. Since no solution has been found, we have removed the symmetry condition 62 = 63. Again 
we do not find any non zero eigenvector E = — 4/3A^ 1 5 E whose singular coefficients satisfy \at\\ — \a%\ = 
|a 3 |. 

6. Conclusion 

The above numerical experiments show that Hypothesis II. II is certainly wrong. The first consequence 
is that we still don't know the shape of the singular set of a minimizer in the situation of Figures 1,2. 
One possibility is that the true singular set is topologically equivalent to the one of Figure 2 but with 
edges tangent to 7 at 7(1) (see Figure 9 below). 

For the moment, this new hypothesis is a conjecture. If it were true, the singular set of a blow-up 
limit at 7(1) would be a half plane and one may expect that the associated global minimizer would be a 
cracktip (type (v)). In this case there is no need to add a new type of global minimizers to the existing 
list in order to explain Figure 9. However one may wonder if there exists a global minimizer whose 
singular set is locally diffeomorphic to the one of Figure 9. Such a global minimizer would not be blow-up 
invariant. 
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c 

Figure 8. Surface: values of |a§, {5 h )\/\a\(8 h )\ for h w 1/40. Black line: values of 
K EF (S h )\/K EF (S h )\ for h = 1/250. 



Another consequence of this negative result is that taking blow-up limit at 7(1), we cannot discrimate 
a surface with a smooth boundary and the surface above. Thus it seems now more difficult to use the 
information on global minimizers to deduce some regularity for the singular sets of minimizers. 




Figure 9. 
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